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Abstract 

We assess the theoretical uncertainties on the inclusive J /if/ production cross section in the Color Evaporation 
Model (CEM) using values for the charm quark mass, renormalization and factorization scales obtained from a fit to 
the charm production data. We use our new results to provide improved baseline comparison calculations at RHIC 
and the LHC. We also study cold matter effects on J '/if/ production at leading relative to next-to-leading order in the 
CEM within this approach. 
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Because the charm quark mass is finite, the total charm production cross section can be calculated in perturbative 
QCD. However, there are large uncertainties due to the choice of quark mass, factorization scale and renormalization 
scale [ 1 1. Typical lower limits of the factorization and renormalization scales are half the chosen charm quark mass 
|Q~||2). In this case, the factorization scale is below the minimum scale of the parton densities. In addition, for 
renormalization scales below 1 GeV, the strong coupling constant a s becomes large and the perturbative expansion 
is unlikely to converge. Thus we seek a set of physically defensible mass and scale parameters that reduce the cross 
section uncertainty. Because the J /if/ cross sections are calculated with the same set of mass and scale parameters as 
open charm production in the Color Evaporation Model (4), we also place limits on the J /if/ cross section calculated 
in the CEM for the first time. See Ref. |3 ] for full details. 

The charm quark mass we employ in our calculations is the Particle Data Group (PDG) value based on lattice 
determinations of the charm quark mass in the MS scheme at /i = m: mini) = 1.27 ± 0.09 GeV |5|. We fit the 
factorization and renormalization scale parameters to a subset of the fixed target total charm production data with 
250 < £beam ^ 920 GeV. The data were evaluated and adjusted to the values we employ in our fits using the most 
up-to-date branching ratios for the measured decay channels in Ref. [6]. We also include data from both PHENIX Q 
and STAR ||8j |9) at -^s NN = 200 GeV. We neglect unknown next-order uncertainties which could be large for charm 
where the mass is relatively small and 0(a).) corrections could be significant. 

The best fit, including the STAR data presented at Quark Matter 201 1 |8 |, yields the parameter values m c = 1.27 
GeV, /J-F/m = 2.1^^ and fXn/m = 1.6^ [j. We show the ;f 2 /dof fit contours on the left-hand side of Fig. [T] for 
Ay 2 /dof = 0.3, 1 and 2.3. The one standard deviation uncertainty in the fitted value of iip/m (fi R /m) was taken as the 
maximum extent of the Ay 2 /dof = 1 contour along the iif/m (yu«/m) axis. The one standard deviation uncertainty in 
the total cross section is the range of cross sections resulting from all combinations of fif/m and fin/m contained within 



Using the final STAR data point |9| in the fitting changes the upper and lower limits on /if/m by 8% and 4% respectively, while the limits on 
/ln/m change by less than 1%. 
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the A^- 2 /dof = 2.3 contour. The A^ 2 /dof = 0.3 contour is to guide the eye. Note the narrow range in pn/m relative to 
the much broader p F /m range. The uncertainty on p F /m is larger and very asymmetric. There is a greater uncertainty 
on the upper limit than the lower limit because there is a much greater change in xg(x,p 2 F ) at lower factorization scales 
than when p F » p.Q, the minimum scale of the parton densities. 
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Figure 1 : (Left) The ^ 2 /dof contours for a fit to the fixed target data and the PHENIX and STAR 201 1 cross sections at V* = 200 GeV. The best 
fit values are given for the furthest extent of the &x 2 = 1 contours. (Center left) The energy dependence of the charm total cross section compared 
to data. The central value of the fit is given by the solid red curve while the dashed magenta curves and dot-dashed cyan curves show the extent 
of the uncertainty bands, corresponding to the furthest extent of the A^ 2 = 1 contours, see text for details. The solid blue curves in the range 
19.4 < ^ < 200 GeV represent the uncertainty obtained from the extent of the A^ 2 = 2.3 contour. Our calculations are compared with the ALICE 
inclusive single muon data from heavy flavor decays 1 1 1| at the LHC for pp collisions at s/s = 7 TeV in the center right and right panels. (Center 
right) The contributions to the pj distributions in (a) divided into rapidity bins, from top to bottom: 2.5 < y < 2.8 (solid red); 2.8 < y < 3.1 (solid 
blue); 3.1 < y < 3.4 (dashed red); 3.4 < y < 3.7 (dashed blue); and 3.7 < y < 4 (dot-dashed red). The top curves are shown at their calculated 
value, the others are scaled down by successive factors of 10 to separate them. (Right) The sum of contributions to the rapidity distribution are 
compared with the FONLL parameter set for charm centered at m = 1.5 GeV (solid red) and our results with m = 1.27 GeV (dashed black). 

The center left panel of Fig. [T] shows the energy dependence of the total charm cross section for the fits with the 
corresponding uncertainty based on results using the one standard deviation uncertainties on the quark mass and scale 
parameters. If the central, upper and lower limits of fi^p/m are denoted as C, H, and L respectively, then the seven 
sets corresponding to the scale uncertainty are {(ji F /m,p F /m)} = {(C, C), (H, H), (L, L), (C, L), (L, C), (C, H), (H, C)}. 
The upper and lower limits on the PDG value of the charm quark mass are 1.36 and 1.18 GeV. The uncertainty band 
can be obtained for the best fit sets by adding the uncertainties from the mass and scale variations in quadrature. The 
envelope containing the resulting curves, 



f max — Ccent ('T'/i.max 'T'cent) 2 (<7"m,max CTcent) 2 ? 

C"min — Ccent — ^ i.^~fi,mm ^"cent) 2 (*7~m,min Ccent) 2 > 

defines the uncertainty. The uncertainty bands are shown for two cases: the fiducial region delineated above, similar 
to Ref. [2|, and including the most extreme cases (pi F /;«,//«/;«) = (H,L) and (L,H). The difference between the 
outer magenta curves, which include these extremes, and the cyan curves, which do not, is very small. Therefore, it is 
reasonable to neglect the extremes. We also show the result for a one standard deviation uncertainty in the total cross 
section obtained from the A^ 2 = 2.3 contour in the blue lines. We have also added the 2.76 and 7 TeV total cross 
sections obtained by the ALICE collaboration in pp collisions 1 10|, not included in our fits. The calculations are in 
rather good agreement with the data. 

We use the FONLL approach 1 2 1 to calculate the heavy flavor semileptonic decay kinematic distributions to com- 
pare to single lepton spectra which include B decays as well as D decays. The B — > fi and B — > D — > p bands 
are calculated with the same fiducial set of parameters as in Ref. The D — > p band is calculated for our best fit 
parameter set. The center right and right panels of Fig. [Tjcompares our calculations with the ALICE single muon data 



(1) 

(2) 
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in the forward rapidity region, 2.5 < y < 4 JTT). The data are given for 2 < pj < 12 GeV and separated into five 
rapidity bins, each 0.3 units wide, as shown in the center right panel of Fig.[T] The calculations agree well with the 
measurements over the entire pj range. On the right-hand side of Fig. [T] we present our results in the dashed curves 
as a function of rapidity, integrated over the same p T range as the data, 2 < p T < 10 GeV. We also show the rapidity 
distribution obtained using the FONLL charm parameter set with a central charm quark mass of m = 1 .5 GeV in 
red. The /^-integrated ALICE data agree well with both calculations. The results with the fitted charm parameter set 
narrow the uncertainty band without sacrificing consistency with the measured data. 

We now turn to a treatment of quarkonium production within this same framework. In the CEM, the quarkonium 
production cross section is some fraction, F c , of all QQ pairs below the HH threshold where H is the lowest mass 
heavy-flavor hadron. We fit F c to the forward (integrated over x F > 0) J/if/ cross section data on only p, Be, Li, C, 
and Si targets. In this way, we avoid uncertainties due to ignoring any cold nuclear matter effects which are on the 
order of a few percent in light targets. We also restricted ourselves to the forward cross sections only. 




Figure 2: (Left) The uncertainty band on the forward J lip cross section calculated based on the cc fit. The dashed magenta curves and dot-dashed 
cyan curves show the extent of the corresponding uncertainty bands. The dashed curves outline the most extreme limits of the band. (Center left) 
The J/if/ rapidity distributions compared to data from 7 TeV [12] (red points and band) and 2.76 TeV [13] (blue points and band). (Center right) 
The forward pj distributions (2.5 < y < 4) are also shown. No additional scaling factor has been applied. A (JcL) kick of 1.49 GeV 2 (7 TeV) and 
1.41 GeV 2 (2.76 TeV) is applied to the pj distributions. (Right) The scale variation of RdAu with the central EPS09 set (blue) compared to the 
EPS09 variation for the central parameter set (red). 

We use the same values of the central charm quark mass and scale parameters as we found for open charm to obtain 
the normalization Fq for (m,pffm,pn/m) = (1.27 GeV, 2.1, 1.6)). We determine Fc only for the central parameter 
set and scale all the other calculations by the same value of F c to obtain the extent of the J /if/ uncertainty band. The 
result is shown on the left-hand side of Fig. [2] 

The ALICE 2.76 and 7 TeV inclusive J /if/ rapidity and forward pj distributions (2.5 < y < 4) are shown in the 
center left and center right plots of Fig. [2] The rapidity distribution at y[s — 7 TeV is flat over several units of rapidity. 
The calculated rapidity distribution at 2.76 TeV is not as broad and the agreement with the data is rather good although 
the midrapidity point remains high relative to the central value of the calculation. The agreement of the calculated 
p T distributions with the forward rapidity data is quite good with the exception of the lowest p T points where the 
calculated distributions turn over more quickly than the data. 

Finally, in the right-hand panel of Fig. [2] we show the ratio Raau for 200 GeV d+Au collisions at RHIC to NLO in 
the total cross section with EPS09 NLO shadowing. The red band shows the variation with respect to the 3 1 EPS09 
sets while the blue band indicates the mass and scale variation around the central set. We note that, with the best fit 
set, the scale uncertainty is now less than the uncertainty due to shadowing. 

We have narrowed the uncertainty band on the open heavy flavor cross section and, in so doing, have also provided 
a realistic uncertainty band on J/if/ production in the Color Evaporation Model. The central result, m = 1.27 GeV, 
p.f fm = 2.1 and pr/iii = 1.6, is quite compatible with previous calculations using a 'by-eye' fit to the data with 
m = 1.2 GeV, p F /m = p R /m = 2 ll4lf!4ll. 
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While the fits have been made by comparing the calculated NLO charm production cross section to available data 
at fixed-target energies and at RHIC, they are in good agreement with the extracted total charm cross sections at the 
LHC. The same parameter set also provides good agreement with the distributions of single leptons from semileptonic 
heavy flavor decays at RHIC and the LHC. The limit on the width of the uncertainty band is now set by the uncertainty 
due to bottom quark production and decay. 

We have used the same fit parameters in the calculation of J/tfr production in the color evaporation model and 
have thus provided the first uncertainty band on J/tfr production in this approach. The energy dependence of the total 
J /if/ cross section that results is a good match to the data up to collider energies. The pj distributions are also in 
good agreement with the data from RHIC and the LHC. In future work, we will use our new parameter set to place 
limits on the contribution of B meson decays to J/ iff production and will also study cold nuclear matter effects on J/tfr 
production in more detail at next-to-leading order. 
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